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The quadratic B-spline curve fitting for the shape of log cross sections 
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Abstract: This paper describes a new method for simulation of the cross section shape of log. The self-developed MQK3102 log shape 
recognizing machine was used to acquire the finite discrete sampling points on the cross section of log and those points were fitted with the 
quadratic B-spline parametric curve. This method can clearly stimulate the real shape of the log cross section and is characterized by limited 
sampling points and high speed computing. The computed result of the previous curve does not affect the next one, which may avoid the 
graphic distortion caused by the accumulative error. The method can be used to simulate the whole body shape of log approximately by 
sampling the cross sections along the length direction of log, thus providing a reference model for optimum saw cutting of log. 
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Introduction 

The contour characteristic of log is extremely complex. To 
recognize the log shape by traditional dimensional measuring 
method may produce errors inevitably between theoretical saw 
cutting and the actual cutting (Tong et al. 1997). With further 
studies on the model of log shape as well as the extensive use of 
computer graphics processing technology in recent years, it has 
been possible to simulate accurately the outline of log on 
computer (Gong 2001). These theories and methods are vitally 
significant to the saw cutting optimization (Ma et al. 2005) and 
the directional cutting of fibre (Ma 2005). The key for simulating 
the outline of log is to establish a proper mathematic model to 
make the fitted graph conform farthest to the actual contour 
characteristics of log. In this study, we adopted the quadratic 
B-spline parametric curve to fit the cross section shape of log, so 
as to provide software and model for the numerical control 
processing of directional cutting of fibre. 

Acquisition of sampling points on log cross section 

We acquired the sampling points on the cross section of log by 
using the self-developed MQK3102 log shape recognition 
machine. Eight detecting pendulum bars are arranged in circle on 
the detecting frame at the center section of the machine. A roller 
with very short diameter is installed at the bottom end of each 
pendulum bar to contact the surface of log, and an angle 
transducer is installed on the mandrel at the fixed end of each 
pendulum bar. When the log moves to the detecting frame, one 
of its ends pushes the eight pendulum bars, and then the eight 
contact rollers fall onto the log. The fluctuation on the surface of 
log causes the pendulum bars to swing, thus leading to the 
change of the output values from the angle sensor. 


Foundation item: The research is supported by Project of National Natural 
Science Foundation of China (30571455) and National “948” Project 
(2005-4-62) 

Biography: REN Hong’e (1962-), Female, Professor, Information and 
Computer Engineering College of Northeast Forestry University, Harbin 
150040, P. R. China. Email: renhonge 163@ 163 .com 
Received date: 2005-12-24; Accepted date: 2006-01-21 
Responsible editor: Chai Ruihai 


Fig. 1 and Fig. 2 show the abridged general view of the 
rectangular detecting frame. The O is an origin of world 
coordinate system in Figs..The horizontal direction along its 
bottom margin is the axis X, and the vertical direction along its 

left margin is the axis Y. Here we suppose that D x and D are 

respectively the horizontal ordinate and the vertical ordinate of 

the pendulum bar’s fixed end. d represents the projection of 
the pendulum bar to the detecting frame plane XOY. CC is the 
intersect angle between the pendulum bar and the negative 
direction of axis X. Therefore, the coordinate value of the 
sampling points on the cross section of log can be defined by 
Equation (1) as follows: 

x - D r +dcosa 

< ( 1 ) 
y -D y — dsincc 

Suppose that the length of each pendulum bar is / , and its 

tilt angle is0, then the projection can be obtained from the 
Equation (2): 

d — l cos 6 (2) 

In Equation (1) and Equation (2), D x , D , CC and / are 

determined by the designing variables of MQK3102 log shape 
recognition machine, and is acquired real-timely by each 
angular transducer. 

Establishment of mathematic model 

B-spline curve is the improved Bezier curve. As to the Bezier 
curve, changing one of the control points will affect the shape of 
curve completely, which is not in favor of the selective 
modification to curve. Moreover, in most cases the curve 
approximated by Bezier can not reflect the features of the 
characteristic polygon (Sun 2000; Sun et al. 1995). 

The parametric equation of quadratic B-spline curve (Sun 
2000) is: 







REN Hong’e et al. 


151 



The expanded form is expressed with matrix as follows: 
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After expanding equation (5), the parametric equation of 
quadratic B-spline curve is obtained according to the descending 
powers of parameter t. 


1 1 2 1 

x(t) = (-x 0 -x l +-x 2 )t + (x, - X 0 )t + - (x 0 + X, ) 

1 1 2 1 

y(0 = (-y o -y\ +-y 2 ) ( +(y i -yo) t+ -(yo +y i) 



From Equation (4) we can get: 



Fig.l Front view of detecting frame 



Fig.2 Side view of detecting frame 


The three vectors P Q , P { and P 2 are resolved respectively 

into P,(x, , JA) which is in X and Y directions on 

two-dimensional space, then the component form is expressed as 
Equation (5): 


/>(<>) = 2 (/>„+/>,) 

P(P = \{Px+Pi) 

Then differentiate the function P{t) to first order: 


P'(t)= X -(2t 1 0) 
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And from equation (8) we can get the results below: 



p'm = p ] -p« 

p'a)=p 2 -p t 


Therefore, we can conclude from Eqs. (7) and (9) that the start 
point and the end point of curve P is the midpoint of P$P X and 

P\P 2 respectively. Moreover, curve P is tangent with P$P X 

and P\P 2 , just like what is shown by Fig.3. 
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Fig.3 Quadratic B-spline curve determined by three points 


From what is mentioned above, we can deduce that to the 
quadratic B-spline curve controlled by Yl (n > 3) points, every 
two neighboring curve segments are tangent at their junction 
point, which satisfies the conditions of both P f (1) = P i+l (0) 
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and p, (i) = p, +1 (0) at one time. Thus, the quadratic 
B-spline curve has fine continuity and smoothness. 

Realizing the Curve Fitting for the Outline of Cross 
Section of Log 

According to equation (6), the parametric equation of the z'th 
quadratic B-spline curve is obtained. It is arranged as follows: 


V(0 = [(-V-1 ~ x i + \ x <+i ) 1 + ( x ; “ x i- 1 )] 1 + y (*«-i + x i ) 

y t ( 0 =[(y y,--i ~y> + 2 y ‘ +l P + U ~ y i- 1 )] 1 +\ (J'.-i + yt ) 


/ - 1 , 2 , • • • 


( 10 ) 


Supposed that the eight sampling points acquired from 
MQK3102 log shape recognition machine are Pq,P\ 

respectively, which are used as the control points during the 
process of quadratic B-spline curve fitting. The parameter 

t (0 < t < 1) is divided into 100 equidistant intervals, that is 

At=0.01. After that we can work out the positions of the fitting 
points according to equation (10) when t is valued differently. 
And then join these fitting points end to end, so the smooth 

segment C- is decided with the three sampling 

points p t _j p, p M . When l is valued from 1 to 6, we can 

obtain six smooth segments of quadratic B-spline curve, and 
every two neighboring smooth segments of quadratic B-spline 
curve are tangent at the junction. 

In addition, because of the shape of the cross section of log is 

the closed, we add two other control points p% and p 9 , and also 
make them meet the conditions of p% — p Q and 
p 9 — p x .That is, we can determine the smooth segments 

C 7 and C 8 which decided by P^PiPo an d PiPoP\ 

respectively according to the method mentioned above. By now, 
we have realized the curve fitting for the outline of cross section 
of log. Here we take one of the cross sections of the log for 
example. The coordinate of eight sampling points along the 
circumference are (290, 365), (210, 301), (200, 222), (270, 159), 
(347, 160), (400, 192), (432, 267) and (408, 326) (The unit is 
millimeter). So after processing in the computer, the fitting 
results are showed by Fig.4. 

Conclusion 

We use the self-developed MQK3102 log shape recognizing 
machine to acquire the finite discrete sampling points on the 
cross section of log and fit those points by the quadratic B-spline 
curve. Such a method can clearly stimulate the real shape of the 
log cross section, does not need many sampling points, and the 
computing speed is high. 

The computed result of the previous curve does not affect the 


next one, which may avoid the graphic distortion caused by the 
accumulative error (Tong et al. 1996). 

As is shown by Fig. 4 (c), if we collects certain cross sections 
in the length direction on the log, and follow the way mentioned 
above, the three-dimensional outline of log may be simulated 
approximately. The closer the sample interval is, the better the 
simulating effect will be. 




c 


Fig. 4 Fitting results 

a. The shape of cross section before curve fitting; b. The shape of cross 

section after curve fitting; C. The effect of simulation after getting certain 

cross sections in the length direction on the log 
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